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?h ! Abstract 

The transverse momentum Qt distribution is computed for inclusive Higgs boson production 
' at y/S = 1.96 TeV. We include all-orders resummation of large logarithms associated with emis- 

sion of soft gluons at small Qt- We provide results for Higgs boson and Z* masses from Mz 

> . 

to 200 GeV. The relatively hard transverse momentum distribution for Higgs boson production 
CN ' suggests possibilities for improvement of the signal to background ratio. 
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I. INTRODUCTION 



Elucidation of the dynamics responsible for the breaking of electroweak symmetry is one 
of the primary goals in particle physics during this decade. In the current run-II of the 
Fermilab Tevatron, the upgraded Collider Detector Facility (CDF) and DO detectors will 
search for the neutral Higgs boson (h), the vehicle of symmetry breaking in the standard 
model (SM), as well as the Higgs bosons of the minimal supersymmetric standard model 
(MSSM). The task is challenging. The production cross sections and branching fractions into 
channels favorable for detection are relatively small, and the backgrounds from competing 
processes are large 111. 

nn 

The fully inclusive gluon-gluon fusion subprocess |2|, |3J gg — > hX supplies the largest cross 
section at Tevatron energies. Its rate is a factor four or so greater than that of the next most 
important subprocess, associated production qq — > hVX, where V = W, Z . In the region of 
modest Higgs boson masses, < 135 GeV, the combination of gluon fusion production with 
the dominant decay process h —>■ bb should be attractive, but its promise is compromised 
by the overwhelming production of bb pairs from strong interactions background processes. 
For m/j > 135 GeV, the channel h — > WW* becomes dominant, and gluon fusion followed 
by h — > WW* provides a potential discovery mode at the Tevatron. 

Calculations of the expected differential cross sections for production of the signal and 
backgrounds are important for evaluation of the required measurement accuracies and detec- 
tor performance. Estimations of the expected transverse transverse momentum distributions 
can suggest selections in this variable that should improve background rejection. In this pa- 
per, we concentrate on gg — ► hX and discuss the behavior of the Higgs boson transverse 
momentum distribution in the region of small and intermediate values of Qt- We expect 
our predictions to be directly pertinent in the region rrih > 135 GeV, but we also entertain 
optimism that experimenters may find clever ways to make use of the dominant channel 
gg — > hX in the region rrih < 135 GeV. 

In the gluon fusion process, production occurs through triangle loops of colored 
(s)particles that couple to the Higgs boson and to gluons. In the SM, the most relevant 
contribution is from a loop of top (t) quarks. The coupling of gluons is simplified in the 
limit of large top quark mass m t |4J . The m t — > oo approximation is valid to an accuracy of 
~ 5% for rrih < 2m t jSj]. Within this approach, the total cross section for gg — > hX, is known 
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to next-to-next-to-leading order (NNLO) accuracy 0,Q|. The large m t approximation serves 
well also for the transverse momentum distribution when rrih < 2m t and the Higgs boson 
transverse momentum Qt is less than m t jq. The next-to-leading order contributions are 



computed in Refs. 

When Qt is comparable to rrih, there is only one hard momentum scale in the perturbative 
expansion of the cross section as a function of the strong coupling a s , and fixed-order 
computations in perturbative QCD are expected to be applicable. However, in the region 
Qt <S m h, where the cross section is greatest, the coefficients of the expansion in a s depend 
functionally on logarithms of the ratio of the two quantities, rrih and Qt, [oc 8 /n) ln 2 (m|/Q^). 
The relevant expansion parameter in the perturbative series is close to 1, and straightforward 
fixed-order perturbation theory is inapplicable for Qt <C rrih. All-orders resummation is the 
established method for mastering the large logarithmic coefficients of the expansion in a s and 
for obtaining well-behaved cross sections at intermediate and small Qt 3, Q Q| • Using 
renormalization group techniques, Collins, Soper, and Sterman (CSS) devised a b space 
resummation formalism that resums all logarithmic terms as singular as (1/Q T ) ln n (m\/Q T ) 



i used widely lor comput, 

HQ HQ 0,0,0 



for computations of the transverse 
and other processes. 



when Qt — * 0. This formalism has been used wide 
momentum distributions of Higgs bosons 

In this paper, we use the b space resummation approach discussed in some detail in our 
Ref. [2^ to derive predictions for the Qt spectrum of Higgs boson production at the Fermilab 
Tevatron energy VS = 1.96 TeV. We resum the large logarithmic terms associated with soft 
gluon emission to all orders in a s obtaining well defined predictions for the full range of Qt- 
We employ expressions for the parton-level hard-scattering functions valid through first- 
order in a s , including contributions from the glue-glue, quark-glue, and quark- ant iquark 
incident partonic subprocesses. We present differential cross sections for masses of the Higgs 
boson that span the range of present interest in the SM, from rrih = Mz = 91.187 GeV to 
rrih = 200 GeV. To illustrate interesting differences, we also provide results for the Z boson. 
In this paper, we also show the integrated Qt distributions cr(<5r > QTmin) for rapidity 
within the interval \y\ < 1.0. This distribution indicates immediately what fraction of the 
Higgs boson cross section may be sacrificed if a selection is made on Qt for background 
rejection. 

Our predictions are presented and discussed in Sec. II. Conclusions and a discussion of 
background rejection are found in Sec. III. 
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II. PREDICTIONS 

Our predictions are based on the formalism described in detail in our Ref. |20[. For 
the perturbative A g and B g functions in the Sudakov factor, we make an expansion valid 
through second order (n = 2) in the strong coupling a s . For the short-distance coefficient 
functions C a ^b used to compute the modified parton densities, we retain the expansion 
through n — 1, i.e., to first order in a s . We remark that the functions A g , B g , and C a ^ 
have well behaved perturbative expansions, free of large logarithmic terms. We use a next-to- 
leading order form for a s (fi) and next-to-leading order normal parton densities <p(x,fi) 

b^- 
In the fixed-order perturbative expressions that enter the "V" function, we select a fixed 



renormalization/factorization scale /i = K^m\ + Q\, with k = 0.5. For the resummed term, 
we take as our central value /i = c/b, with c = 2e _7B ~ 1. We examine sensitivity to the 
choice of /i by computing cross sections with other choices of this scale. The extrapolation 
into the non-pert urbative region of large b is accomplished with the form devised by Qiu 
and Zhang 23}. 

In Fig. [TJ we show the predicted differential Qt distributions. We present results for 
three choices of mass of the Higgs boson, = Mz = 91.187 GeV, = 125 GeV, and 
m/j = 150 GeV (where the WW* decay channel is dominant). In all cases, the solid lines 
represent the predictions of resummation at next-to-leading order (NLO) accuracy and, 
for comparison, the dashed lines show the results of resummation at leading-order (LO) 
accuracy j^ . 

Two points are evident in the comparison of Z boson and Higgs boson production, with 
m/j = Mz, Fig. ^c) and (d). The peak in the Qt distribution occurs at a smaller value of 
Qt for Z production. At y — 0, the curve peaks at Qt ~ 3.0 GeV for Z production and 
at Qt ~ 7.3 GeV for Higgs boson production. Second, the distribution is narrower for Z 
production, falling to half its maximum by Qt ~ 9.0 GeV, whereas the half-maximum for 
Higgs production is not reached until Qt ~ 21 GeV. The physics behind these important 
differences is that the larger QCD color factors produce more gluonic showering in the 
glue-glue scattering subprocess that dominates inclusive Higgs boson production than in 
the fermionic subprocesses relevant for Z production. After all-orders resummation, the 
enhanced showering suppresses the large-6 (small Qt) region more effectively for Higgs boson 
production. 
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FIG. 1: Differential cross sections for Higgs boson production at V^S* = 1.96 TeV at (a) Q = mh = 
125 GeV, (b) Q = nih = 150 GeV, and (c) Q = Mz- The NLO prediction is shown as a solid line 
and the LO prediction as a dashed line. In (d) we show analogous results for Z boson production. 
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For Higgs boson production, there is a substantial difference in the predictions at NLO 
and at LO, but the differences are slight for Z production. The role of the C function is 
res pon sible for these effects. The first two coefficients in the expansion of C are published [17I 
[rsl. llo| - In powers of (a s /7r), the C^> coefficients for the gg and qq cases are 



Cfl g {z) = 5{l-z)\ 



\c F {l-z). (1) 



In the qq case relevant for Z production, there is a strong cancellation, 2ir 2 — 24, reducing 
the effect of the NLO contribution. For Higgs boson production, the analogous term is 
additive, 2tc 2 + 11. The constant terms —24 and 11, respectively, are associated with virtual 
corrections at NLO and can be traced to the interference between the tree and one-loop 
diagrams. The 2n 2 — 24 term is effectively the finite piece of the one-loop correction to the 
qq —>■ Z vertex while the 27r 2 + ll term results from the finite piece of the one- loop correction 
to the effective gg — > h vertex. 

Integrating our Qt distributions over all Qt and rapidity y, we obtain the total cross 
sections as a function of Higgs boson mass shown in Fig. |3 We may compare these cross 
sections to the values listed in Table 3 of Ravindran et al y\ based on their fully inclusive 
calculation, without resummation. Their results are obtained with different parton densities 
and are quoted at \/~S = 2 TeV, making the comparison somewhat imprecise. The integral of 
our NLO resummed cross section lies between the NLO and NNLO inclusive results, slightly 
above the NLO values, for all masses shown, as would be expected qualitatively. Our 
calculation includes the ingredients and Bjp of a NLO calculation with, in addition, 
other ingredients such as of the NNLO total cross section. 

Comparison of the predicted Qt distributions for Higgs boson production at different 
masses of the Higgs boson shows that the peak of the distribution shifts to greater Qt as 
rrih grows and that the distribution broadens somewhat. At y = 0, the peaks are centered at 
about Q T = 7.3, 8.2, 8.7, and 9.5 GeV for m h = M z , 125, 150, and 200 GeV, respectively. 
The change of the Qt distribution with irth can be examined quantitatively with plots of the 
mean value < Qt > and of the root-mean-square < Q\ >2. We find that < Qt > grows 
from about 19 GeV at = Mz to about 28 GeV at rrih = 200 GeV. The curve may be 
approximated with a straight line over this range, with < Qt >— 0.079m/j + 12 GeV. The 
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FIG. 2: The integral of the resummed cross section over all Qt and rapidity for Higgs boson 
production at t/S = 1.96 TeV as a function of Higgs boson mass at NLO (solid line) and LO 
(dashed line) accuracy. 

root-mean-square of the distribution grows from about 27 GeV to about 38 GeV, reflecting 
the broadening of the Qt distribution with m^. 

For comparison with Higgs boson production, we quote our predictions for Z production: 
< Q T >= 13 GeV and < Q\ > 1 ' 2 = 20 GeV. The difference < Q T > - < Qf >~ 6 GeV at 
rrih = Mz is a manifestation of more significant gluonic radiation in Higgs boson production. 

In anticipation of a discussion of selections on Qt to improve signal to background ratios, 
it is instructive to examine production of Z* bosons, with mass greater than Mz- For Z* 
masses above Mz we assume that the dominant qq production model is unchanged except for 
the difference in mass of the Z*. We find that the change in the Qt distribution with mass 
is much less significant than is seen for the Higgs boson. For example, the peak position at 
y = increases to only 3.4 GeV at Mz* = 200 GeV from its value of 3 GeV at Mz- 

We introduce a new distribution for Higgs boson production that represents the integral 
of the differential cross section for all Qt greater than a minimum value, and for rapidity 
integrated over the interval \y\ < 1.0. The results are shown in Fig. El At QTmin — 10 GeV, 
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FIG. 3: (a) Cross sections for Higgs boson production at %/£ = 1.96 TeV integrated over the ranges 
Qt > QTmin and \y\ < 1-0 for three values of the mass of the Higgs boson: mh = 125 GeV, 150 
GeV, and mh = Mz- In (b), we show analogous results for Z boson production and, for ease of 
comparison, the curve for Higgs boson production at = Mz rescaled so that the cross section 
at QTmin = is the same as that for the Z case. 

the integrated distribution has dropped to 40% of its total for Z production but to only 
65%, 71%, and 74% for Higgs boson production at mh = M z , 125 GeV, and 150 GeV, 
respectively. The irreducible background for Higgs boson decay to pairs of W's arise from 
the QCD annihilation subprocess qq — > WW. This subprocess has the same initial state 
structure as Z(Z*) boson production. The harder Qt spectrum for the glue-glue dominated 
Higgs boson production shows that the signal to background ratio can be enhanced if Higgs 
bosons are selected with large Qt- 

Choices of parameters are made in obtaining our results, including the renormaliza- 
tion/factorization scale fi and the non-perturbative input. Our default value fi = c/b with 
c = 2e~ lE ~ 1 provides a scale that varies inversely with the impact parameter b. This 
selection has the virtue that logarithmic dependence on u is removed from the coefficient 

nn 

functions C [13, 20]. The integration in b space is dominated by the region near the peak in 
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the b distribution, b ~ 0.08 for = 125 GeV. Using the conjugate relationship Qt ~ 1/b, 
we note that the typical hard scale fi is therefore < fi >~ 12.5 GeV (~ O.lm/j). To examine 
sensitivity to the selection of /i, we consider other choices that are independent of b. Taking 
a value as large as fi = 2m^ produces changes in the magnitude of da/dydQj- at the peak 
position of no more than about 20%, much less than the difference between the NLO and 
LO results in Fig. Uncertainties associated with non-perturbative physics in the region 
of large b are at most 1 to 2% depending on the size of the power corrections j^. 



III. CONCLUSIONS 



Discovery of the Higgs particle is essential to shed light on the mechanism of electroweak 
symmetry breaking. The partonic subprocess g+g — > hX dominates Higgs boson production 
in hadronic reactions when the Higgs boson mass is in the expected range < 200 GeV. 
The two-scale nature of the production dynamics, with mass and transverse momentum 
Qt both potentially large, and the fact that the fixed-order perturbative QCD contributions 
are singular as Qt — > 0, necessitates all-orders resummation of large logarithmic contribu- 
tions in order to obtain meaningful predictions for the Qt distribution particularly in the 
region of modest Qt where the cross section is largest. We perform this resummation of 
multiple soft-gluon emission using an impact parameter 6-space formalism. The resummed 
Qt distributions are determined primarily by the perturbatively calculated 6-space distri- 
butions at small b, with negligible contributions from the non-perturbative region of large 
b. 

In this paper, we present predictions for the Qt distributions of Higgs boson and Z 
production at v^S" = 1.96 TeV. At the same mass, = M z , the predicted mean value 
< Qt > is about 6 GeV larger for Higgs boson production than for Z boson production. 
For the Higgs boson, < Qt > grows from about 19 GeV at rrih = Mz to about 28 GeV at 
vrih = 200 GeV, and the root-mean-square < Q T > 1//2 from about 27 GeV to about 38 GeV. 

Searches for the Higgs boson in its decay to WW require an excellent understanding of 
the production characteristics of both the signal and backgrounds. The irreducible back- 
grounds arise from the QCD annihilation subprocess qq — > WW. Of interest to us is the 
expected Qt dependence of the ratio of signal to irreducible background. In this paper we 
provide predictions for the transverse momentum distribution of the signal. The annihila- 
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tion subprocess qq — > WW has the same initial state structure as Z(Z*) boson production. 
As we show, the Qt spectrum of Z* production is predicted to be significantly softer than 
that for Higgs boson production. We suggest therefore that a selection of events with large 
Q^ w should help in improving the signal to background ratio. The placement of the cut 
requires, of course, appropriate optimization to maintain signal significance. In the mass 
range rrih < 135 GeV where h — > bb is the leading decay mode, the QCD subprocess gg — > bb 
supplies the dominant background. Since the signal and background are both produced 
from gg scattering, we expect them to have very similar Qt dependences at = m b i, and 
means other than selections on Qt will be necessary for background suppression in this mass 
interval. 
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